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We study the critical behavior of the d = 3 Ising model with bond randomness through extensive 
Monte Carlo simulations and finite-size scaling techniques. Our results indicate that the critical 
behavior of the random-bond model is governed by the same universality class as the site- and 
bond-diluted models, clearly distinct from that of the pure model, thus providing a complete set of 
universality in disordered systems. 
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Understanding the role of impurities on the nature of 
phase transitions is of great importance, both from exper- 
imental and theoretical perspectives. First-order phase 
transitions are known to be significantly softened under 
the presence of quenched randomness [l|-lH| i whereas con- 
tinuous transitions may have their exponents altered un- 
der random fields or random bonds 0, 0]- There are 
some very useful phenomenological arguments and some, 
perturbative in nature, theoretical results, pertaining to 
the occurrence and nature of phase transitions under the 
presence of quenched randomness [1, H, ISl] ■ Historically, 
the most celebrated criterion is that suggested by Har- 
ris . This criterion relates directly the persistence, un- 
der random bonds, of the non random behavior to the 
specific heat exponent ap of the pure system. According 
to this criterion, if Up > 0, then disorder will be relevant, 
i.e., under the effect of the disorder, the system will reach 
a new critical behavior. Otherwise, if ap < 0, disorder is 
irrelevant and the critical behavior will not change. Pure 
systems with a zero specific heat exponent (ap = 0) are 
marginal cases of the Harris criterion and their study, 
upon the introduction of disorder, has been of particular 
interest (lo| . The paradigmatic model of the marginal 
case is, of course, the general random two-dimensional 
(2D) Ising model and this model has been extensively 
debated [uj. 

Respectively, the 3D Ising model with quenched ran- 
domness - which is a clear case in terms of the Harris 
criterion having a positive specific heat exponent in its 
pure version - has also been extensively studied using 
Monte Carlo (MC) simulations and field theo- 

retical renormalization group approaches [20l - f2^ . Es- 
pecially, the diluted model can be treated in the low- 
dilution regime by anal ytic al perturbative renormaliza- 
tion group methods [23l425| . where a new fixed point in- 
dependent of the dilution has been found, yet for the 
strong dilution regime only MC results remain valid. His- 
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torically, the first numerical studies of the model sug- 
gested a continuous variation of the critical exponents 
along the critical line, but it became clear, after the 
work of Ref. [l3| , that the concentration-dependent crit- 
ical exponents found in MC simulations are the effec- 
tive ones characterizing the approach to the asymptotic 
regime. Note, here, that a crucial problem of the new 
critical exponents obtained in these studies is that the 
ratios P/v and ^/v occurring in finite-size scaling (ESS) 
analysis are almost identical for the disordered and pure 
models. In fact, for the pure 3D Ising model, accu- 
rate values are [26]: = 0.6304(13), (3/iy = 0.517(3), 
j/v = 1.966(3), and a — 0.1103(1). Respectively, for the 
site- and bond-diluted model, the most accurate sets of 
asymptotic exponents {jS / v,^ / v,v) have been given by 
the extensive numerical works of Ballesteros et al. fly| 
and Berche et al. [Hi are: (0.519(3), 1.963(5), 0.6837(53)) 
and (0.515(5), 1.97(2), 0.68(2)). 

The above estimates of critical exponents provided ev- 
idence that the 3D Ising model with quenched uncorre- 
lated disorder belongs to a single universality class, dis- 
tinct from that of the pure model, as also indicated by 
the Harris criterion, independent of the considered disor- 
der distribution. Yet, in a very recent paper Murtazaev 
and Babaev using MC simulations and ESS meth- 
ods on the site-diluted model the above view was contra- 
dicted and these authors suggested that the model has 
two regimes of critical behavior universality, depending 
on the nonmagnetic impurity concentration. Motivated 
by the above contradictions and the great theoretical in- 
terest of the existence of universality classes in disordered 
models, we have chosen to investigate the 3D Ising model 
with bond disorder in order to compare all these three 
kinds of disorder (site-, bond-dilution and bond disorder) 
and to verify whether these lead to the same set of new 
critical exponents, as would be, in principle, expected by 
universality arguments . In this contribution we show 
that, indeed, the 3D Ising model with quenched, uncor- 
related bond disorder belongs to the same universality 
class as the site- and bond-diluted models, defining in 
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this way a complete universality class in disordered spin 
models. 

In the following we consider the 3D bond-disorder Ising 
model whose Hamiltonian with uncorrelated quenched 
random interactions can be written as 



(1) 



where the spin variables Si take on the values — 1,+1, 
(ij) indicates summation over all nearest-neighbor pairs 
of sites, and the ferromagnetic interactions Jy > follow 
a bimodal distribution of the form 

'P{J^,) = I W^J - Jl) + 5{.h, - J2)], (2) 

where Ji -I- J2 = 2, Ji > J2, and r — J2/J1 reflects the 
strength of the bond randomness. Additionally, we fix 
2kB/{Ji + J2) = 1 to set the temperature scale. The 
value of the disorder strength considered throughout this 
work is r = 1/3. 

Resorting to large scale MC simulations is often nec- 
essary, especially for the study of the critical behavior of 
disordered systems. It is also well known that for such 
complex systems traditional methods become inefficient 
and thus in the last few years several sophisticated al- 
gorithms, some of them are based on entropic iterative 
schemes, have been proven to be very effective [1^. The 
present numerical study has been carried out by apply- 
ing our recent and efficient entropic scheme [1^. In this 
approach we follow a two-stage strategy of a restricted 
entropic sampling .30] based on the Wang-Landau (WL) 
algorithm [3l|. As we do not wish to reproduce here the 
details of our implementation, we give only a brief discus- 
sion on the nature of the WL method. The usual WL re- 
cursion proceeds by modifying the density of states G{E) 
according to the rule G{E) fG{E) and initially one 
chooses G{E) = 1 and f = fo = s. Once the accumula- 
tive energy histogram is sufficiently flat, the modification 
factor / is redeflned as: fj+i = y/Jj, with j — 0, 1, . . . 
and the energy histogram reset to zero until / is very 
close to unity (i.e. / = e^""' « 1.000 000 01). Once / is 
close enough to unity, systematic deviations become neg- 
ligible. However, the WL recursion violates the detailed 
balance from the early stages of the process and care is 
necessary in setting up a proper protocol of the recursion. 
In spite of the fact that the WL method has produced 
very accurate results in several models, it is fair to say 
that there is not a safe way to access possible systematic 
deviations in the general case. This has been pointed out 
and critiqued in a recent review by Janke [32|. However, 
from our experience and especially from our recent stud- 
ies on 2D and 3D disordered spin models , the WL im- 
plementation followed in these papers has produced ex- 
cellent results, enabling the discrimination between com- 
peting theoretical predictions on that model. Using this 
combined approach we performed extensive simulations 
for several lattice sizes L e {8,16,24,32,40,64}, over 
large ensembles of random realizations of the order of 
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T =4.4219(49) 
V = 0.6843(67) 
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FIG. 1: Shift behavior of several pseudocritical temperatures 
as defined in the text. The error bars reflect the sample-to- 
sample fluctuations. The inset shows the FSS of the sample- 
to-sample fluctuations of the pseudocritical temperature of 
the magnetic susceptibility. 



1000 — 3000. Each disorder realization was simulated at 
least 10 — 20 times with different initial conditions to 
improve accuracy. 

It is well known that, extensive disorder averaging is 
necessary when studying random systems, where usually 
broad distributions are expected leading to a strong vi- 
olation of self-averaging [l3, H^- A measure from the 
scaling theory of disordered systems, whose limiting be- 
havior is directly related to the issue of self-averaging 17 1 
may be defined with the help of the relative variance 
of the sample-to-sample fluctuations of any relevant sin- 
gular extensive thermodynamic property Z as follows: 
Rz = {[Z'^]av-[Z]l^)/[Z]l^. Closely related to the above 
issue of self-averaging in disordered systems is the man- 
ner of averaging over the disorder. This non-trivial man- 
ner may be performed in two distinct ways when identi- 
fying the flnite-size anomalies. The flrst way corresponds 
to the average over disorder realizations ([. . .]av) and then 
taking the maxima ([. . .]*„), or taking the maxima in each 
individual realization first, and then taking the average 
([. . .*]av)- Closing this brief outline, let us comment on 
the statistical errors of our numerical data. The statis- 
tical errors of our WL scheme on the observed average 
behavior were found to be of small magnitude (of the 
order of the symbol sizes) and thus are neglected in the 
figures. On the other hand for the case [. . .*]av the error 
bars shown reflect the samplc-to-sample fluctuations. 

We briefly present here our numerical results for the 
3D random-bond Ising model. Figure [1] illustrates in the 
main panel the shift behavior of 6 pseudocritical tem- 
peratures estimated via the second way of averaging dis- 
cussed above, i.e. by taking the average over the individ- 
ual pseudocritical temperatures. The error bars reflect 
the sample-to-sample fluctuations. The pseudocritical 
temperatures considered correspond to the peaks of the 
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y/v= 1.967(3) 
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FIG. 2: FSS of the maxima of the disorder-averaged magnetic 
susceptibility in a log-log scale. The inset shows the limiting 
behavior of the ratio R[x']av ■ 



following six quantities: specific heat C, magnetic suscep- 
tibility X, derivative of the absolute order parameter with 
respect to inverse temperature K = 1/T: d{\M\) /OK ^ 
{\M\H) - {\M\){H) [si] and logarithmic derivatives of 
the first {n — 1), second {n — 2), and fourth (n = 4) 
powers of the order parameter with respect to inverse 
temperature d\n{M'^)/dK = {M'^H)I{M'^) - (H) [H. 
Fitting our data for the whole lattice range to the ex- 
pected power-law behavior [TJlau = Tc + bL~^/^ , where 
Z stands for the different thermodynamic quantities men- 
tioned above, we find the critical temperature to be 
Tc = 4.4219(49) well below the critical temperature of 
the pure model and the estimate v — 0.6843(67) for the 
critical exponent of the correlation length. This value is 
in excellent agreement with the values 0.6837(53) and 
0.68(2) given by Ballesteros et al. tl(| and Berche et 
al. 

Using the above sample-to-sample fiuctuations of the 
pseudocritical temperatures and the theory of FSS in 
disordered systems as introduced by Aharony and Har- 
ris [sj] and Wiseman and Domany [l^, one may fur- 
ther examine the nature of the fixed point controlling 
the critical behavior of the disordered system. Accord- 
ing to the theoretical predictions [13, [13 ' the pseudocrit- 
ical temperatures of the disordered system are dis- 
tributed with a width d[T^]av, that scales, in the case 
of a new random fixed point, with the system size as 
S{[Tz]av) ~ L-^/". In the inset of Fig. [U we plot these 
sample-to-sample fiuctuations of the pseudocritical tem- 
perature of the magnetic susceptibility. The solid line 
shows a very good power-law fitting giving the value 
0.688(15) for the exponent ly, which is also in very good 
agreement with the value 0.6843(67) obtained via the 
classical shift behavior and the most accurate estimates 
in the literature [H, [l^ . 

In Figs. [5] and |3] we provide estimates for the mag- 
netic exponent ratios of the model. In particular, in 



FIG. 3: FSS of the disorder-averaged order parameter at the 
estimated critical temperature in a log-log scale. The inset 
illustrates the FSS of the maxima of the disorder- averaged 
inverse-temperature derivative of the absolute order parame- 
ter. 



Fig.[2]we present the FSS of the maxima of the disorder- 
averaged magnetic susceptibility [%]*„ in a double log- 
arithmic scale. The solid line presents a linear fitting 
using the total lattice range spectrum, giving the esti- 
mate 1.967(3) for the ratio j/v {[x]*ay i^/"")- The 
inset shows the ratio Rz, where Z = [x*]a-u as a function 
of the inverse lattice size indicating a strong violation 
of self-averaging of the magnetic properties of the 3D 
Ising model with bond disorder. Respectively, in Fig. [3] 
we plot the disorder-averaged magnetization at the esti- 
mated critical temperature, as a function of the lattice 
size L, also in a log-log scale. The solid line is a lin- 
ear fitting {[M]av{T = Tc) ~ L-'^^") giving the value 
/3/i' = 0.516(11). Additional estimate for the ratio /3/i/ 
can be obtained from the FSS of the derivative of the 
absolute order parameter with respect to inverse tem- 
perature which is expected to scale as with the 
system size 135]. Thus, in the corresponding inset of 
Fig. Owe plot the data for d{\M\) /dK averaged over 
disorder as a function of L, also in a double logarithmic 
scale. The solid line is a linear fitting that gives an es- 
timate of 0.518(9) for the ratio /S/i/. Thus, overall the 
values for the ratios of the magnetic exponents are in 
excellent agreement with the best estimates for the site- 
iP/i' = 0.519(3) and 7/1/ = 1.963(5)) [3 and bond- 
diluted {13/iy = 0.515(5) and -//i^ = 1.97(2)) ^ cases 
reinforcing the scenario of a single distinctive universal- 
ity in the 3D Ising model with quenched uncorrelated 
disorder, independent of the disorder distribution. 

Using our estimates for the critical exponents and the 
Rushbrooke relation {a/iy + 2/3/i' + j/i' = Ijv) we es- 
timate the ratio ajv to be —0.085(10). In Fig. [4] we 
present the FSS of the maxima of the disorder-averaged 
specific heat data as a function of the linear size L. The 
solid line shows a typical power-law fitting attempt of the 



4 



1 , 1 , 1 , 

a/v--0.08(l) 


1 ' 1 ' 1 ' 


0,004 




/ 0,003 








p p^- 0,002 








/ 0.001 




\ 




/ 0.000 

, / 1 , 1 , 1 






0.00 

1 


0.05 0.1(1 0.15 0.20 0.25 
1 /L 

1,1, 



I U. I I I I I I I I I I I I I 

10 20 30 40 50 60 70 

L 



FIG. 4: FSS of the maxima of the disorder-averaged spe- 
cific heat. The inset shows the limiting behavior of the ratio 

form [C]*„ ~ L"/"^, that gives a value ajv = -0.08(1). 
This is a further rehabihty test of our numerical method 
and extended simulations. The corresponding inset of 
Fig. m presents the ratio i?[c*]ai. ^ function of the in- 
verse linear size. As in the inset of Fig. [2l also here a 



violation of self- averaging is observed, yet this is, in ab- 
solute numbers, much smaller than that of the magnetic 
susceptibility. 

Summarizing, we have presented in this brief report 
concrete evidence that the critical behavior of the 3D 
Ising model with quenched uncorrected disorder is con- 
trolled by a new random fixed point, independent of the 
way randomness is implemented in the system. This re- 
sult has been obtained through extensive numerical sim- 
ulations and classical finite-size scaling techniques. Par- 
ticular interest was paid to the sample-to-sample fluc- 
tuations of the pseudocritical temperatures model and 
their scaling behavior that was used as a successful al- 
ternative approach to criticality that verified the above 
scenario. Although we acknowledge that Ballesteros et 
al. [l^ and Berche et al. were the first to support 
numerically the above view, we believe that the present 
contribution puts a further significant step on this in- 
tensively debated issue of universality in disordered spin 
models. Very interesting would be also to study more 
complicated models, where disorder couples to the order 
parameter, and one such prominent candidate is the 3D 
random- field Ising model [SG*]. For these type of mod- 
els, the existence of universality classes has been severely 
questioned fSTj . 
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